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ABSTRACT 

detem1n2£? St ? dy T a PP bcation of radio interferometry to sateUite orbit 

t de ? loped a ““Nation tool for assessing interferometric tracking 
SSS?;- Ini ? aUy J u V ie l t001 was applied to the problem of determining optimal interferometric 
station siting for orbit determination of the Tracking and Data Relay Satellite (TORS') 

Subsequently, the Orbit Determination Accuracy Estimator (ODAE) was expanded to model the 

‘‘Nation ■*»)«-. of the 

^termination System (GTDS) with measurement types including not only group and nhase delav 

rate, rngular mtLSSS^SS SSR! * 
satellite measurements. The user of ODAE specifies the statistical properties of error sources 
including inherent observable imprecision, atmospheric delays, station bcation uncertainty and 
measurement biases. Upon Monte Carlo simulation of the orbit determination process ODAE 
calculates the statistical properties of the error in the satellite stale vector mSZ&tiS muameters 
for which a solution was obtained in the orbit determination. V Pieters 


n ^ pres f nts f es ^ ( te from 0DAE application to two different problems: 

W ■™" mnatl ®P of optimal geometry for interferometric tracking of TDRS, and (2) expected 

for u Global Positioning System (GPS) tracking of low-eanhorbit 
t f Ult ^!' C ®? c l ua ®5 # about 0 P timal ground station locations for TORS orbit 
remniM a r ^y radl ?.f terfe 1 i ;? met ry are presented, and the feasibility of GPS-based tracking for 
IRIDIUM, a LEO mobile satellite communications (MOBILSATCOM) system, is demonstrated 
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INTRODUCTION 

A* nart of its effort to assess cost and performance benefits of various emerging technologies, 
fponS a Series ofstudies on the application of radio interferometry to 
^l^acMn^ThoSgh astronomers had used radio interferometry for decades P n °r> itwas not 
until the late 1960s that interferometry was proposed for use in satellite orfelt ,^j™ naU ‘ 

experiment devised by Irwin Shapiro, Alan WWtney,^d o re, veiy ^ satellite in 

interferometric fVLBD measurements were made on the TACSAT I communications , sateime ui 
^synchronous orbit (GEO), and the semi-major axis of the orbit was measured wtth accuracy 
on^dieorder of several hundred meters [1], Subsequent experiments If by 

Jim Rav Curt Knight, and others to determine the posiuon of the Tracking and Data Relay 
WlUte (Ssiyielded accuracy on the order of 75 meters [2]. Such orbit determination 
?ShS«m the extremely high precision of the group delay and phase delay 
observables, make radio interferometry an attractive option for satellite tracking. 

Operational considerations are also a benefit of radio interferometry in satellite orbit 

detenmnation, because the group and phase dday me^urem^ 

Whereas the existing Bilateral Ranging Transponder System (BRTS) is taxing on 1 uixa 
c^ScSions^ Surces, radio intorferome£y can derive its measurements from any signal 
including the signal intended for the TORS user community 

ftp-termination svstem for TDRS would eliminate traffic for tracking on the TDRS transponder. 
Because an interferometric tracking system would be passive, it would place no design constraints 
on the space segment, and it would therefore provide backward compatibility with all generations 
of TDRS. Thus, NASA found radio interferometry to be an attractive technology to pursue 
future TDRS tracking applications. 

NASA sponsored a series of studies to investigate whether an operational radio interferometry 
system couldprOTide TORSorbit determination Srvices (1 ) at lower cost, (2) at greater ^curacy, 
S ^^considerably smaller baselines than BRTS. Contributors to these studies included 
“Irferome^Inc wKere a Small Business Initiative Research (SBIR) contract™, th NASA 
wS executed to demonstrate hardware and software that would provide group delay 
measurements on TDRS with VLBI. CSC performed an assessmentforthe Godd^dSpace 
Plight Center (GSFC) on a variety of TDRS tracking alternatives, including VLBI and Connected 
Flement Interferometry (CEI) systems. The Jet Propulsion Laboratory (JPL) sponsored a senes 
of^wrim^^to^remine CEI accuracy from its Goldstone facility. For its part o the effort, 
MITRE assessed optimal site locations and expected life-cycle costs of an operational 
interferometric TDRS orbit determination system. 

For accuracy assessment purposes, MITRE developed a Monte Carlo simulation tool, the 
OrbU Destination Accuracy Estimator (ODAE), that initially modeled error sources in orbit 
detenrtination with VLBI aiJ CEI systems. In ODAE, the user can specify a sateUite °*^any 
2t of ground stations between which group or phase delay measurements are to be madeandthe 
statistical properties of the error in those measurements. Upon each iteration of the Mo ^Carlo 
simulation theorbit of the satellite is determined based on measurements with errors added, and 
the errors in the resulting satellite ephemerides are recorded. Thus, the user may study the 
statistical properties of the error in the batch orbit determination process resulting from the use of 
r»°™. The initial application of ODAE was to study the effects of 

varvine satellite and measuring station geometries on orbit determination accuracy and to propose 
5H by ™»i° interferometry. The results of that study are presented 

in this paper. 

Subsequent studies for the United States Air Force (USAF) on Space SurveUlance b^twork 
fSSNl accuracy and Global Positioning System (GPS) accuracy led to the expansion of ODAE to 
rare, azimuth andlleUm angle, and satellrte-tosatcUrte measurement 


386 


types. The potential application of GPS to satellite tracking has been under consideration for a 

nmh^ r ° f f? 81 ® &&’ . For . 9 F 9 ^telUtes, which have higher orbits than GPS satellites the 

However fo^FJRO^tSi^ 116 “J* 1 weak signal strength present limitations [4]. 

However, for LEO satellites, which have lower orbits than GPS satellites, visibility and signal 

SiJif 6 greaUy . im P roved ; Attention has been focused recently on LEO satellites in tearena 

llt^ mm r CaUOnS - A , T nber ofLEO MOBILSATCOM systems arecurrentWin 
plannmg or development, including Motorola’s IRIDIUM, Loral’s Globalstar, TRW’s Odyssey 

hfccTw 1 TSTKW?^ T JSS aj0in 5 ve " ture of Unship Partners with William Gates and Craig ’ 

^<RiStSiSSiSSS^' accuracy of GPS tracki " 8 for 


THE ODAE MODEL 

^ batch maximum likelihood orbit determination process applied in the 
* iUOn S ? s *y ( ° TDS) [«]• specifies a reference tree 

fte tnnrs at which measurements are to be made. Given a set of obse^Uons ™ the 

measurements, group or phase delay measurements, or pseudorange measurements) 
ODAE determines the set of parameters (e.g„ state vector, clock offsets, or Sspheri c ’ 

(WAESi-K^r^f 51 f r ** obsen ' all0[ ' s - Upon each iteration of its Monte Carlo Emulation, 
ODAE injects errors of user-specified statistical properties into various parts of the orbit 

de y™J all0 “ Process. ODAE computes the error of the measured parameters at each iteration 
and at the end of the simulation, ODAE computes the statistical characteristics of the eiror. ’ 

E™ ^ can be modeled by ODAE include inherent measurement imprecision 

atmos P heric delays, “» d clock offsets. The user mustspecify the 
statistical properties of the error sources. Trajectory propagation schemes for dynamic orbit 

9°^ f J om ^e two-body approximation to numerical integration of the fully 
foSf T h" S 0f A detailed mathematical specification of 8 he coordinate frame 

force models, and numerical integration techniques used in ODAE are given in Reference 7 The 

sSr ffiSf *T** from , ^ G P S a PP roach to «** determination is the use ^ of Buhrsch- 
Stoer rational function extrapolation for numerical integration [8, 9]. For the numerical 

“S 1 ^ equations of satellite motion, the Bulirsch-Stoer technique has been shown to 

provide the same precision as more traditional techniques, such as predictor-corrector integration 

ri mtegratl0n ’ but 5- reduced computational cost [7, 10]. For short-teim dynamic 
orbit determination accuracy studies to assess the relative effects of changes in station eeometrv or 

10 app,y ***** " on schemes A°e 

• ...9 DAE . wa f implemented in Mathematica to allow maximum flexibility of the model Since its 
TDR^MFTRF^h 10 ^ P r oWem of optimal ground station siting for interferometric tracking of 
TDRS, MPTOElm applied ODAE to a variety of problems. Most recently, MITRE has proposed 
the use of ODAE for assessment of initial orbit determination accuracy with the HAVE STARE 
a P pb PT 0ns include the assessment of orbit determinafon accuracy for tiilSce 
fl ^ Pr0Ve, ? ent * 0gram ( SSNn> ) for various classes of orbite, L to the 

bvSrMarh^S TT Cy Van ° US SCen , an0S - ^ough computational time is increased 
Mathematica, it allows for very natural representation of the equations of motion 

mencal mtegration schemes, and the batch orbit determination algorithm Also because ODAE 

f Cld L adap “ We *°. a "** ofproble m "Ld^nl°Hifc 
0f 0DAE 10 TORS -0 *» GPS navigatio/for a 
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GROUP DELAY AND PHASE DELAY MEASUREMENT FUNCTIONS 

Consider an interferometric mbit deteimination scenario in which O is the origin 10 J _an Mr*- 
centered inertial (ECI) coordinate system, r is the position vector of a satellite with resjtfct to , 
h and h are the Dosition vectors of two ground stations from which measurements are to be 
mate ^ vectors of the satellite with respect to those ground scions 

as pictured in Figure 1. Theposition vectors r, b t , b 2 , d^ and d 2 are all functions of time. .The 

sum ofariation |osition vector, b„ , and the sateLite position vector ' So^ateTof 

d. is simply the satellite position vector r; therefore, d* = r- b* . If the propagation ra^, c, 01 
tte signal through the atmosphere is known, then the transit time, T k , of the signal from the 
satellite at point P to ground station number k at point B k will be given by 

T k =i|d*| = iV( r " b *)-( r_b *) 


The true group delay, T, between stations 2 and 1 is the differential transit time of the signal 
between these two sites: 

r=T 2 -T l = £(jd 2 | - (d t J) = ^(r-b 2 )(r-b 2 ) - ^/( r - b i)- (r - b t )] 


P 



Figure 1. Illustration of the Interferometric 
Measurement Scenario 


A subtletv of equation (1) is that the satellite position vector, r, the station 1 position vector, 
. ^dAesrationTposition vector, b 2 , are all referenced to different times. If the measured 
signaltmanates fronnhe satellite at time t, then it will arrive at station 1 at time t + T, and it will 


388 


amve at stouon 2 at time t +T+ x, where T is the signal transit time from the satellite to station 1, 
and tis the true group delay between stations 1 and 2. If the satellite position vector r is 
measured at time t, then the station 1 position vector is measured at time t + T, and the station 2 
position vector is measured at time t+T+x. Thus, we write r = r(r), b t = b.(r + T) and 
t> 2 - b 2 (r + T+ x). The group delay equation (1) is, therefore, more properly written as follows- 


T = i Vl r (0 - b 2 {t + T + x)\ • [r(0 - b 2 (r + T + t)] 

~ i Vt r W - b i(' + r )H r (0 -bj(/+ r)] 


( 2 ) 


io»;t, Il a 0 ? AE v. a J 1Ser s P e ?‘. fie * a scenario that includes epoch time, satellite state vector at epoch, 

alut J ld ® of interferometer sites, and times at which measurements are to be 
taken. ODAE must then calculate the true group delay observables from equation (2), but as can 
be seen, the right-hand side of (2) is a function of r. Therefore, ODAE solves equation (2) for x 
iteratively, as described in Reference 7. During the Monte Carlo simulation, ODAE computes 
measured group delay by adding measurement or atmospheric fluctuation errors to the true group 
delay as computed above. The solution of the orbit determination problem on each iteration of the 
simulation, as described in Reference 7, follows the GTDS maximum likelihood estimation 
approach one step of which is the computation of the Jacobian, or matrix of partial derivatives of 
equation (2) with respect to the state vector parameters at epoch. 

• ,L^ 0r P* 1356 delay measurements, ODAE converts phase delay into equivalent group delay. If v 
is the reference center frequency of the phase delay measurement, AT is an integer number of 
signal cycles, and ^ is the true phase delay, then the equivalent true group delay x can be 
computed as follows: 1 


_ Qj¥2nN 
2ttv 


( 3 ) 


This computation can be accomplished so long as the cycle ambiguity N can be determined from a 
p non information about the satellite’s position vector. ODAE allows the user to model clock 
offsets or local oscillator offsets for group or phase delay measurements, respectively. In such 
cases, the offset is taken as an additional parameter in the orbit determination process. 

APPLICATION TO TDRS 

In thissection, the level of GEO satellite orbit determination accuracy that can be attained with 
radio mterferometry is demonstrated, and conclusions about optimal station-satellite geometry are 

!rmoek are applied to recommend optimal ground station siting for orbit determination 

ot IDRS by radio interferometry. 

Radio interferometry with baselines the size of BRTS’s, which are intercontinental, would 
translate tte; high level of observable group delay accuracy into greatly improved TDRS tracking 
accuracy. However, it was NASA’s desire instead to accept only a modest improvement in 
accuracy whfle reducing system cost and ameliorating other operational considerations by greatly 
shortening the baselines. This led naturally to the study of a CEI-based system, where barlines 
are very short Because of the requirement for a CEI system to have interferometer sites 
connected by fiberoptic cable in a temperature-controlled environment, the cost of lengthening 
baselines is very high. We constrained our baselines to 20 km maximum length for the purposes 
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We used ODAE to assess position determination accuracy for a baseline scenario and to 
determine the effects of varying the relative satellite to ground station geometry. Because the 
effect only of relative geometry was to be studied initially, it was not necessary to 

TDRS ephemerides or true potential ground s ^tion /^cations. nJovide three 

geosynchronous with a 4° inclination and a subsatellite longitude of 18 . P r’d 

fndependent baselines across which phase delay could be measured, we constmned ou 
sites to lie on the vertices of a square with a 20 km baseline, as shown in Ague 2. The** 
latitudes longitudes, and altitudes for this reference scenario are given m Table 1. ODAE 
modeled’ simultaneous phase delay measurements acro^ the barlines £j m ^uon 2 to ’ 

station 3 to station 1, and station 4 to station 1 (denoted 2-1, 3-1, and 4-1, respectively), inese 
baselines are illustrated in bold in Figure 2. 


i 

0°W 


i 



Table 1. CEI Station Locations for 
Reference Scenario 


Station 

Number 

Geodetic 
Latitude (°V 

Longitude 
(°E) 

Altitude 

flan) 

1 

45.00000 

0.0000 

0.1 

2 

45.00000 

-0.2545 

0.1 

3 

45.17997 

0.0000 

0.1 

4 

45.17997 

-0.2545 

0.1 


An extension of Alan Whitney’s work [11] shows that the theoretically achievable precision of 
the phase delay observable, Cfy , is given by 


1 

” 2 n{SNR)v 


(4) 


where vis the center frequency, in Hz, sampled by the interferometer. Since the TDRS downlink 
toWhite Sands is centered at 14 GHz, the theoretically achievable precisionof tiiephase delay 
observable is 0.23 picosec. While no TDRS tracking experiments were performed with JPL s 
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CEI equipment at Goldstone, observations were made on natural radio sources at 8.4 GHz to 
assess the precision of the phase delay observable [12]. The statistical phase error, expressed in 
radians, is roughly 1/SNR, and JPL’s experiments at Goldstone demonstrated a typical phase 
error of 0.005 cycles [13]. The achieved SNR was therefore l/(2*x0.005) = 32. At 8.4 GHz, 
relationship (4) predicts a phase delay observable precision of 0.59 picosec. JPL demonstrated 
the standard deviation in the phase delay observable to be approximately 1 picosec [12], which is 
a factor of 1.7 larger than the theoretically achievable value. Applying this factor to the 
theoretically achievable phase delay precision for TDRS, we estimated the practically achievable 
precision to be 0.23x1.7 >=* 0.4 picosec. We took this measurement error to be independently 
normally distributed across each baseline. For the initial study, it was assumed that there were no 
equipment biases, that there were no atmospheric delay errors, that all station were connected by 
fiberoptic cable to one clock and frequency standard, that there were no local oscillator offsets 
between the four stations, and that station positions were known with perfect accuracy. Thus, the 
pure effect of measurement geometry and observable precision on orbit determination could be 
assessed. 

ODAE Monte Carlo simulation of the orbit determination scenario described above with 200 
iterations showed a lcr root-mean-squared (RMS) position vector accuracy of 3.2 km. We also 
assessed the accuracy that can be attained with the use of other combinations of baselines. It is 
practical to have one site in common for all three measurements so that the common site can act as 
the correlation center at which the phase delay observables are generated. For the particular 
satellite and ground station locations in this scenario, selection of three measurements where one 
station is common to each pair (i.e., 2-1, 3-1, 4-1; or 1-2, 3-2, 4-2; or 1-3, 2-3, 4-3; or 1-4, 2-4, 
3-4) results in a lcr RMS position vector accuracy of 3.2 km. Thus, there is no geometrically- 
preferred common site for the measurements. 

The orbit determination scenario described above was the starting point for the assessment of 
the effects of varying interferometric measurement geometry on orbit determination accuracy. 
Since only relative geometry matters, and since it would have been more cumbersome to vary the 
positions of four ground stations, we instead varied the satellite’s initial position vector. First, we 
studied the effect of relative interferometer baseline size on orbit determination accuracy. Satellite 
range from station 1 was varied while keeping the elevation angle and azimuth angle from that 
station constant. Because the baseline sizes are small relative to the range to GEO, the range, 
elevation angle, and azimuth angle from each of the other three stations is close to that of the first 
For this particular orbit determination scenario, range from each site to the satellite is 
approximately 37,850 km, the elevation angle is approximately 39°, and the azimuth angle is 
approximately 155°. As shown in Figure 3, the smaller the range to the satellite for a constant 
baseline length (or, equivalently, the longer the baselines across which phase delay is measured 
relative to the range to the satellite), the greater the position vector accuracy. 

Next, we assessed the effect of satellite azimuth angle on orbit determination accuracy. The 
azimuth angle of the satellite at station 1 in the original scenario was varied while keeping the 
range and elevation angle from that station constant Figure 4 shows the variation in position 
determination accuracy with satellite azimuth angle. The results indicate that for a configuration of 
four mterferometric ground stations at the vertices of a square, position error is maximized when 
the satellite s azimuth angle is an integer multiple of 90°, and position error is minimized when the 
satellite s azimuth angle is an integer multiple of 45°. 

Finally, we assessed the effect of satellite elevation angle on orbit determination accuracy in 
this scenario. The elevation angle of the satellite at station 1 was varied while keeping the range 
and azimuth angle from that station constant As can be seen in Figure 5, for this particular orbit 
determination scenario, position error increases monotonically with elevation angle. Thus, based 
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Figure 5. Position Error vs. Satellite Elevation Angle 


on the criterion of minimizing ephemeris error due only to error in the phase delav measurement 

<tegene ra ™Ste3te “ ***** ^ “ d ‘ he sceoario becomes 

ratafnSST ra !f throu ^ h ^Posphere show that errors in predicting signal propagation 
«te elevation angle decreases [14, 15]. Moreover, errors in predicting propagation 

CEI svstt^THl P ^enc fluctuauons tend to be the dominant error source in overall accuracy for 
Cbl systems [15]. Thus, we sought to determine the optimal elevation angle for CEI * 

measurements with consideration of both measurement error and tropospheric delay error. 

We modeled tropospheric fluctuations between each interferometer site and the satellite as 
being independently normally distributed. The assumption of independence is based on the fact 
at water vapor cells can be of several kilometers in diameter, and so tropospheric delav errors 
mo^ eaCh S1 ^. can V 1 fac 5 ** independent For an elevation angle of 20° and for a 100 second 
measurement duration, the magnitude of the standard deviation in tropospheric delav error was 

SS2K to l L 4 piC0Sec ;i, 4] - Ovation angle dependence of SS devTationl 

f i U T ^ SqUare r0 ° f t0f ** Structure Unction calculated in 
• A lso ’. unde f ? e assumption of independent errors along each station-to-satellite 
path, the variance in phase delay error for a particular measurement pair will be Ae sum 5the 
variances along each path. Since elevation angles are roughly equal along each path the standard 
delation of the phase delay error will, therefore, be V2 ti mis the standafdSation Sne 

SJrtT A ? Se ^“^P^sand the results in Reference [14], we computed the values of 
standard deviation m phase delay error due to tropospheric fluctuations shown in Table 2. 

n. 1 i 0r f Varying ® levation ^gies. we used ODAE to model error due to tropospheric 

fluctuations as well as inherent phase delay imprecision. The resulting lcr position etrors are 

Flgure 6 - As can be seen, the optimal satellite elevation angle is approximately 30° In 
tite conclusions section of this paper, we show how these results can be applied to optimally 
siting a CEI system for TDRS orbit determination. PP pUaUy 
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Table 2. Standard deviation in phase delay error 
due to tropospheric fluctuations as a function or 
elevation angle 


Elevation 
Angle O 

Tropospheric Delay 
Error (picosec) 

10 

7.5 

20 

5.7 

30 

4.6 

40 

3.9 

50 

3.3 

60 

3.0 



APPLICATION TO GPS 

Studies of GPS accuracy for a precision landing system and other GPS-related smdies for the 
Department of Defense led to the application of ODAE to GPS accuracy problems ODAE models 
the P GPS navigation problem in much the same way as it models the satellite orbit^tormination 
nmhlern In the case of GPS navigation, a “station location” is the position of a GPS satellite at a 
particular time GPS almanac data are loaded into ODAE and propagated to the desired 
S5££S times Unknown receiver position can be on the surface of the earth or on a 
satellite and in the case of the former, the satellite orbit is also propagated to the desired 
measurement time OD AE tests the visibility of satellites in the GPS constellation from the 
receiver position, and if the number of visible satellites exceeds the number of available channels 
in the receiver, ODAE determines the optimal subset of satellites for measurement 
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To demo,B,rate -S-cy. ~ 

T*®# 1 * designed to consist of 66 satellites in orbits with a 
11 S8°^[o1!?L 0f 7 ’- 43 kpj 16 !- Satellites will be divided into six orbital planes spaced bv 
3 .58 in nght ascension of the ascending node. Orbital inclinations will be 86 4° and y 

SSkMf.!?, 3 ' .*?*“. ?■ ta V** equally^ S* pl^a^ijacem 

out of phase with one another. We selected an arbitrary satellite from the 

We^^OlTpLlT < n r hl Jf« nsi< ? oflhe "<xte 31.58“, inidal mean anSy 0° for Sysis 
fWtKn determm€ : number of GPS satellites visible to this IRIDIUM satellite a/a 

function of time, and we used ODAE to determine position accuracy as a function of time* 

^25 I J° St *£5* G ^? almanac data C2/i3/9 4 at this time of this study) included the full 

C.R^ h^am° n .^ 6 satelllte f For the purposes of assessing GPS satellite visibility we assumed a 
Si^ ldth ofa PP roxi mately 27°, which is the angle subtended bytheSUm Ss 
n£ «Lf} < f aUSe the *!f? m wi l lh «**. than 27°, it wo“d S 

would be required to ^W^v^thtTnecemry g^,°^TT^C w^ight MviagTwo^kfb^^ 00 ^ 

Gps to * 

IRtnlfru. 7 ,S?o ,S If num ! Kr of 0PS sateUites visible as a function of time from the reference 
M . Ca " be .«*»t f“- foot hours, fewer than four GPS sateUte m Se 
thC IRIDIUM “*>“* wo»M be 0.83. However, srith 
because only thni visibKuta wo "“^^7 P ° sWoni " g S “ vices ,hrough GPS 


V 



Time (hours) 


Figure 7. GPS Satellites Visible from an IRIDIUM 
Satellite Over a One-Day Period 
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We next computed position accuracy as a function of time for IRIDIUM tracking by GPS. 
Fl J£ of precision (PDOP) across a 24-hour period for mstanmneous 
£ fi^Xr SrpSts where only three GPS satelUtes were visible, we assumed that an 
accurate clock was available and that the navigation 

expected. Figure 8 shows poorest accuracy at the times when fewer GPS stellites are visioie. 
Except for those times, baseline PDOP appears to be on the order of 2 or 3. 

Finally we used ODAE to compute dynamic orbit determination accuracy for IRIDIUM 
traclringby GPS^Figure 9 shows PDOP across a 24-hour period where me^^mente are taken 
on the hour added to the previous pool of measurements, and processed in bate . 

PDOP dee reases below a value of one, and after 10 hours, a value of roughly 0.5 is 
S.TF^ow^c^acv situation where the precision of the pseudorange measurement is on 

mnmJM Wm^£5g|SS r 

of 0.5x(3 meters) = 1 .5 meters. Such accuracy is likely to be sufficient, even with IR1D1U 
stringent formationkeeping requirements [16]. 



Figure 8. Position Accuracy over a One-Day 
Period for an IRIDIUM Satellite 


CONCLUSIONS 

by siting a CEI system in eastern Louisiana or western Mississippi. 
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Figure 9. Dynamic Orbit Determination Accuracy over a 
One-Day Period for an IRIDIUM Satellite 


In the second part of this report, we assessed GPS availability and accuracy for tracking a 
LEO satellite. In particular, for IRIDIUM, where the proposed constellation consists of satellites 
m near-polar orbits at altitudes of 785 km, instantaneous orbit determination accuracy is available 
at the level of 9 meters (lo), and long-term dynamic orbit determination can reduce errors to the 
level of 1.5 meters. Because GPS receiver equipment has the potential of offering reduced weight 
and cost by comparison with traditional TT&C equipment, GPS provides an attractive tracking 
alternative for LEO satellites such as IRIDIUM. 
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